


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1984 


Analysis of the failure mechanism of 
Shipboard structural members exposed to fire. 


Hiddemen, Robin Lee 


Massachusetts Institute of Technology 


http://ndl.handle.net/10945/19135 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 


/ (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist : Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

; | LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


| 
| 








DEPARTMENT OF OCEAN ENGINEERING 


MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
CAMBRIDGE, MASSACHUSETTS 02139 


ANALYSIS OF THE FAILURE MECHANISM OF 
SHIPBOARD STRUCTURAL MFMBFRS EXPOSED TO FIRE 


ish 
ROBIN LEE HIDDEMEN 


OF COURSF XIII-A 
SM(NA & ME) JUNF 1984 






bia HO Lareiee - 
Bia Pay - 
i 


y 
> 


° —_ 


= 





ANALYSIS OF THE FAILURE MECHANISM OF 
SHIPBOARD STRUCTURAL MEMBERS EXPOSED TO FIRE 


by 
ROBIN LEE HT ODEMEN 
} 


R.S. Monmouth College 
ler 73 


Suomitted to the Department oft 
Ucean Engineering 
in Partial Fulfillment of the 
Requirements of the Degrees of 
JCEAN ENGINEER 


and 


MASTER OF SCIENCE IN NAVAL ARCHITECTURE 
AND MARINE ENGINEERING 


at the 
Peo oseaUsei tS INSTITUTE OF TECHNOLOGY 
June 1934 
(C) Massacnusetts Institute of Technology 
1784 
The author hereby grants to the United States Government 


ang its agencies permission to reproduce and to distribute 
copies of this thesis document tn whole or in part 


Ane gi 4 
WS272NS 





ANALYSIS OF THE FAILURE MECHANISM OF 
SHIPBOARD STRUCTURAL MEMBERS EXPOSED TO FIRE 


by 
ROBIN LEE HIDDEMEN 


submitted to the Department af Ocean Engineering on 11 
May 1°89 in partial fulfillment of the requirements of 
the Degrees of Ocean Engineer and Master of Science in 
Naval Architecture and Marine Engineering 


RESTRACT 


The objective of this study is to analyze the behavior of 
a typical ship structural member expased to fire. 


Included in this work are chapters pertaining to _ the 
modelling of the fire process and the structural member toa be 
analyzed. After determining the expected temperature exposure 
Ot a shipboard structural member, a heat transter analrtis 15 
carried out and the resulting load is subsequentiy used in the 
failure analysis. 


iiememasvoc thrust of thie study is the analysis art the 
failure mechanism of the collapse of ant! beam. Two methods 
Of analysis developed by Kecman, and Hayduk and Wierzbicki, 
are summarized and applied to the case of a fire load. The 
recuced moaments due to tire exposure are determined and 
finally recommendations are offered tor the use of this 
information in ship structural design. 


Thesis Supervisor: Tomasz Wierzbicki 
Title: Professor of 
Applied Mechanics 





€2LENOWLEDGEMENTS 


IT would like to acknowledge the invaluable assistance and 
quidance af my advisor, Protessar Tomasz Wierzbicki. 
Additionally I would like acknowledge my husband, Louis T. 


Seaega tor his continued support during my three years at MIT. 





la fe ot eoncemes 


jm tle Paqe 
Abstract 
AcKnowledgements 


met Of Fiqures 
List of Tables 


Pier educt i on 


rs 


‘Chapter 


Model of Structure 
Backqround 
Load Transmission 
Selection of Material 


Chapter 


ho to No bh 


Whe 


Me@ael Of Fire Precess 
Charseteristics of a Shipboard Fire 
Development ot a Time Temperature Curve 


Chapter 


GW W® W 
ij — 


Heat Transter Analysis 
Conduction Heat Transfer 
Selutian 
Results 


Chapter 


BRL 
W tom 


Analysis of the Failure Mechanism 

Kecman”’s Method 

Je-l1.1 Maximum Bending Strength 
Actual Bending Collapse Mechanism 
Theoretical Collapse Mechanism 
Nominal Energy Absorbed During 

Hinge Ratation 

The Final Form of the Hinge Moment 
and Wierzbicki’s Method 

Collapse Modes 
au Isometric Transformation Modes 

2 Hinge Collapse Modes 

2.1.3 Assembled Collapse Modes 


Chapter 


Ol cn 
-— 


Gi Un cn 
Sr 
bw 


ch 
h 
a) es Oe a) 
e ie 
hx 
a - 
~~ c Ol 
ra 


on of Results ofr the Two Methads 


Conclusions and Recommendations 
if Summary 
2 Recommendations 


Chapter 


Table of Centents 


Page 
References Fe 
Appendix A Heat Transter Analysis Program Wag 
Appendix B Kecman’s Analysis Program 101 


Appendix C Hayduk and Wierzbicki’s Analysis Program 104 





&o BL 
hWwn- 


oC OS Gh 


gnguauo vo 


ae ard 
Died o 


Slst ot figure = 


Typical Moment Curvature Relations 


Typical Midships Section 

Typical Shipboard Structural Beam Elements 
Dimensions of Model to be Analyzed 
strength Temperature Relations of 

ASTM-ASS Steel 


Three Periods ot Fire 

Duration of Full Developed Fires 

Average Temperature of Fire Gases 

for a Typical Compartment 

Time Temperature Curves for Various Types 
or Fire 


aimeecon+!Guration 


Thermal Conductivity of Plain Carbon Steel 


Results of the Heat Transfer Analysis 
Dimensions of the Beams Analyzed 


Deformation Pattern 

Paeiacracal Surface TU 

Model of Hinge Collapse Mechanism 
The Anaqle of Rotation alang GK 
Results af Kecman’s Analysis 
Illustrations of Inmextensional and 
Extensional Deformation Paths 
Assembling Folding Modes into More 
Complicated Patterns 

Deformation Patterns 

Model of T Section 

Falding Mechanisms 

Crushing Case 

Extensional Collapse Mode - Definition 
of Geometric Parameters 

Bending Case 

Definition of the Angle of Rotation 
Moment Displacement Relation 
Comparison of Approximate and Exact 
solutions 


Results of Hayduk and Wierzbicki’s Analysis 


Comparison of the Two Methods 


Steel 


Specitic Fire Loading of Shipboard Spaces 


Properties 


List of Tables 


Paqge 


18 
22 


Chapter 1 


INTRODUCTION 


Traditionally, Marine Engineers have handled the effects 
of fire on shipboard structures through the use of damaqe 
control procedures. Ideally, a structure should be designed 
to withstand anticipated fire loads. ici oGbieecti ve of this 
work is to apply the methods of post critical analysis of 
thin-walled structures developed by Mechanical Engineers to 
analyzing fire loads on board ocean engineering structures, 


specifically ship structures. 


A ship structure is a highly redundant structure. If the 
load bearing capacity of one member is diminished, then the 
surrounding structural members carry more of the applied 
load. Shipboard fires are characteristically highly localized 
phenomenon which may effect only a small percentage cof the 
entire structure. By examining the effect of fire on a single 
structural member, one will gain insight into the behavior of 
the entire structure. Ultimately, ship structures may be 


designed to withstand the effects of typical shipboard fires. 


It jis the intent of this study to examine the effect of a 


fire on a typical shipboard structure, namely the I beam. The 
major thrust will be the analysis of the failure mechanism of 
the beam. Two methods will be introduced which analyze the 
post yield behavior of an I = beam. Figure eet shows 
graphically how ae fire load diminishes’ the mechanical 


properties of a beam. 


The material yields upon reaching the maximum moment in 


Figure 1.1. Consider the following relation: 


Stiffness = Moment/Curvature., 


Beyond the yield point, stiffness is less than or equal to 
zero. Normally a structural member remains in the elastic 


range, and the moment remains less than M , the fully plastic 
=) 


yield moment. However, through exposure to fire, the 
stiffness will decrease so rapidly that the moment will § fall 


below a safe range. 


The ultimate objective of this study is to show how much 
the strength of an I beam is’) diminished when exposed to fire. 
To this end, the structure and the fire process will be 
modelled; and the effect of temperature on the mechanical 
praperties of a selected material will be given. After these 
models have been developed, a heat transfer analysis will he 
performed. As aresult of this analysis, it will be shown 


that bending collapse is the primary mechanism ot structural 
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failure. Next, the failure mechanism analysis methods will be 
fully developed. Finally, recommendations will be offered for 
applying the results of this work to the design of an entire 


ship structure. 





Chapter 2 


MODE! OF STRUCTURE 


This chapter will develop the motivation towards choosing 
an I beam as a typical shipboard structural element. 
Additionally, a material will be selected; and the temperature 


dependent properties of that material will be given. 


fet Back oround 


The principal aspects of ship structural design most 
commonly employred by Marine Engineers invalves the applicatian 
of rules established by Societies such as the American Bureau 
of Shipping. These rules are based upon the primary loads 
experienced by a ship during its fii eercanes These loads 
include: still water bending moments and shearing forces; wave 
induced moments and shearing forces; impulse loads; thermal 


loads; and dynamic loads. 


These loads which are unique to ocean gaing vessels are 


mecem ditticult to accurately opredict. Therefore, ship 
structural design is largely a process involving guesses 
derived from intuition or empirical data. The primary stress, 
strain and displacement behavior of the ship structure are due 
to displacements limited to those parallel to the unloaded 
midsurface of the plate. Thus, traditionally hull strength is 
determined through an analysis of the midships section. 


Figure 2.1 illustrates a typical midships section. 


fee toad [ransmission 


Today, Ship Structural Engineers have adopted a made! 
which is motivated By considering the loads transmitted to the 
structure. Loads transmitted to a stiffened plate (by 
slamming, etc.) are in turn transmitted to the stiffeners and 
from them to the point supports, and so on. Eventually a load 
transmitting grillage is derived. Next, boundary conditions 
Will be used to define if the spans are treated as a series of 
isolated points or continuous beams. In most cases the 
boundary points will be those which can be described as simple 


or clamped supports. 


Figure 2.2 shows examples of typical members of a ship’s 


structure. For the purposes of this study, an I beam will be 
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analyzed as its behavior 1s considered representative of the 
open beams in Figure 2.2. Figure 2.3 gives the dimensions of 
the beam which hereafter will be used as the model to be 


analyzed. 


The first moament of area of the I beam will be required 


im the analysis. 





tae 
vs zadA = t.(b - Ea (Es yee ; W (221) 
A = Sectian Area 
t = Thickness of the flanges 
+ 
t = Thickness of the web 
Ww 
a = Width of flanges 
b = Height of web 
iZ 
Assume t is negligible. Hence, the first moment of area may 


be approximated by: 
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2.3 Selection of Matecial 


Several materials have been used throughout shipbuilding 
history starting with wood, and progressing to iron and 


steel. Today several grades of steel and even titanium are 


imam 
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used in the shipbuilding industry. One of the mast  comman 
grades of steels used is ABS Grade A. Material Scientists, 
however, have not tested ABS Grade A in order to determine the 
temperature dependent material properties. Ordinary 
structural steels such as ASTM-A36, have been widely used as a 


test specimen. Table 2.1 compares the properties af the two 


steels. 

Table 274 

Steel Praperties 

ASTM-A3S ABS Grade-A 
Soarben, max “~A El, eS G326 
Phasphoaraus, max % Q.04 0.04 
Sulpher, max ~% 0.05 0.04 
Yield Stress, Kei 36.00 394.00 (mins 

It is considered a canservative assumption that an 

analysis for ASTM-AS6 can be used without sacrificing 


mecuracy, «ref. 1). 


One af the major cantributians ta the loss of strength is 
the temperature dependence af the mechanical properties af the 
materials under consideratian. Several researchers have 
predicted the deteriacration af steel structurai elements at 
elevated temperatures. A search af the Jiterature reveals 
good agreement among thase whom have studied the temperature 


dependence of material properties. 





Harmathy and Stanzak (ref. 2) have carried out several 
experiments in order to provide to Design Engineers, primarily 
to those concerned with fire endurance of structural elements, 
the elevated temperature mechanical properties of common 
commercial steel] ArOoGuGc tc. The temperature dependent 
expression for Young’s Modulus of Elasticity, E, ‘psi? “for 


ASTM-A3SS steel is given here: 


S De 
Be=—230*20- ==. 9ST 


ho 
G) 
Nw 


oO 

T = Temperature of heat expasure, Fahrenheit 
Steel loses strength as it is heated. Figure 2.94 illustrates 
the temperature dependent tensile strength of ASTM-ASS6 steel. 


a a) 
Notice that between 1000 and 1200 F the tensile strength 


falls below the usual design stress. 


For the purposes of this study, linear dependence ahs 
temperature is assumed. The following equation will be used 


in this analysis. 


on = —2@ 5) tae 7,000.0 
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Chapter 3 


MODE! OF FIRE PROCESS 


A survey of the literature reveals several methods of 
characterizing a fire. Harmathy offers a numerical technique 
which is summarized in Section 3.1. Most Fire Engineers use a 
time temperature curve to define the characteristics of a 
fire. The development of an acceptable time temperature curve 


f= Given in Section 3.2. 


Gee ee baracteristics of _ a Shipboard Fire 


The heat exposure on a shipbeard structure resulting from 
Pere is a function of several identifiable amd incidental 
Votameiples, (rem. 3). The three primary factors contributing ta 
the severity of fully developed fires is the flow rate of air 
into a Burning compartment, ‘compartment ventilation); the 
total fire toad; and the thermal characteristics of the 


compartment boundaries. 





Numerous theoretical etudies have determined 


time-temperature curves of typical fires based on these 
factors, ‘ref 3). Even with simplifying assumptions, these 
factors are only calculable With a computer. Harmathy 


introduces three easily calculable fire severity parameters 
which can be assumed to be accurate predictions of the 
Sontributing factors ta a fully developed fire. These 


parameters are defined as follows: 
newer ation of the fully developed fire, Tt; 


2. The overall penetration flux, gq, that is the heat flux 
absorbed By the compartment boundaries, averaged 


spatially over the boundary surfaces, and; 


3. The average temperature of the compartment gases, 


Caverage “fire temperature”), T , averaged spatially over 
g 


the period of full development. 


It is useful at this time to introduce the three periods 
of freely burning compartment fires. See Figure 3.1. The 
focus of this study i¢ on the period of full development, 
since 80-90 percent of the fuel energy is released during this 


ae | Od. 


A compartment fire is a fire which is enclosed in a shell 
of arbitrary shape. Typically» on board a ship, ventilation ts 


secured ta the compartment with the fire, and free 
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communication to the surrounding spaces is not allowed. Fire 
load is a numerical expression which is a function of the 
available fuel in the compartment. The units of specific fire 


load is pounds of wood per square foot, (ref. 1). 


Ta compute the fire load of a particular compartment, the 
products of the weights in pounds of the material, and the 
heat capacity are summed for all the combustible materials in 
the compartment. The resulting expression, which has units of 
BTU’s must be divided by 8000 BTU/Ib. This result will provide 


the same available eneraqy as the equivalent weight of wood. 


Equation €3.1) gives the equivalent weight of wood in 


pounds divided by the deck area of the compartment. 


= (3.13 
F= : WH. /H AL 


Zz 
F = Specific fire load, I]b/tt 


W = Weight of individual materials, Ib 
H = Heat capacity of individual materials, BYU/1b 
H = Heat capacity of wood, BITU/IE 
Ww 
2 
A = Area of deck, ft 
7 


Table 3.1 gives typical fire load values for various 


- 249 - 


shipboard compartments, (ref. 4). 


Tabl eae. i 
specific Fire Loading of Shipboard Spaces 


los./square feet 

Control Spaces 
Wheelhouse/Char troom AS) 
Fire Control Stations Ces) 


Escape Routes 


Corridors 1.3 
Stairway Enclosures 1.0 
Accommodation Spaces 
Staterooms 
Fire Resistant Furnishings oO 
Combustible Furnishings oo 
Public Spaces 
Fire Resistant Furnishinas: 
Lounges, Restaurants, etc. Seu 
Ferry Vessels ao 
Combustible Furnishings 2.0 
Sanitary Spaces 0.4 
Service Spaces 
Galleys Ag, 
Pantries no heating appliances) 4.0 
Food Concession (ferry vessels- 
no combustible stowage) io 
Workshops Pos.0 
Storerooms - Combustibles 10.0 
Storerooms - Cleaning gear anly 3.0 
Ship’s Laundry Ora 
Laundry - Private Use te 
Main Machinery and Cargo Spaces 10.0 
Auxiliary Machinery Spaces Zall 
Tanks and Voids 0.0 


Tt has Been determined By several researchers that the 
flow rate of air into a Burning compartment, (ventilation) is 
just as critical to the severity of the fire as the fire 
load. fhe determination of compartment ventilation depends on 
how the rate of entry of air compares with a critical value. 
The rate of entry of air into a compartment can be expressed 


aS5 


U =c an (3.2) 
a lowow 

ee ow fate of air, |b/hr 

a 

z 

A = Area of window, ft 

W 
h = Height of window, ft 

Ww 

of e 

C = 2450 Iib/Ft =r 


The critical value is defined as: 


Oar 5 CoOFAL C€3.3) 


Ze 


C= 23.) eet “nr 


O= Material constant which depends on the furnishings of 


the compartment. For typical living spaces, 





toes < 8 =< 0.70, 


The air flow factor is expressed as follows: 


le eee (3.49) 
a 


inte, a fire is ventilation controlled if ¢ ¢ 1.03; and a 


fire is fuel surface controlled if € >= 1.0. 


The average rate of Burning of combustibles, R, has been 
determined during the periad of full development based upon 
empirical data. For § >= 1.0 

Ee CEA ee, 
a COPA, 
Za 
Cees 1b tt eahir 
= 
Bepee <« 1.0 
R= Onda (3.0) 


The mass flow rate of gaseous products out of the 


monear tment is: 


t 
aw) 
“J 

} 


moe follows that: 


Q = R(0.9326H, + 0.0686H.) ¢ 


3) 
an 
Ne 


where: 


Q = average rate of heat evolution from combustible 


materials 


H =~* 37°60 BIU/IB, heat of combustion of char 
1 


H =~ 11480 BTU/Ib, heat of combustion af volatile 
2 


Gecomposition products 


68 = the fraction of heat released by the combustion of 
the volatile decomposition products inside the 


compar tment 


6 16 calculated Sased upon an estimation ot the average 


flame length, |: 


: 1/3 (Sao) 
& = C, (OFA) 
1/3 
Be — pe ae ie ie oe 2 
P| 
en 63.10) 
5a 
1/3 17s 
C = 90.366 tho 7 lia 


- 23 - 


Then, 


< 
AE Age He 6 =] CScel 
if 2>4 Bs ln Pye (312) 
; Cc Gs 
where H is the height of the compartment. 


5 


The duration of a fully developed fire, +t, is calculated 


as follows: 


fore t= ¢,/ (3.13) 
2 
C = 0.205 #t Chr)/Ib 
é 
tO cul T= C./06 Co eho 


The other two parameters may be obtained numerically from 


the following equation: 


i ae 4 4 S S5 
=_ —_ —_ —_ _ CS 
q Ae Sag See ed Seas TS )] 

wheres 
A = Total area for the compartment 

t 
Q = Average rate of heat evolution (BFU/hAr) 

a) 

C = Specific heat of compartment gases (BTU/IB- R) 

9 


I 
Nh 
“2 

\ 





Oo 


T = Temperature of combustion gases, R 
Q 
a) 
T += Temperature at time = 0, R 
O 
~8 2 o 64 
pemoeet an—-bOltzman constant, 0.17110 Set Carin Ce 
ma—sempirtiecal factor, 1.05, dimensionless 
= Clo q melo G7 
“gy Wee ele eames elects) ie ie (Seep 
wheres 
fe— emoair ical tactor, 0.7, dimensionless 


| 2 
K = Average thermal diffusivity, ft /hr 


a 
K = Average thermal conductivity, BTU/Ft-hr- R 


Figures 3.2 and 3.3 show the dependence of the averal] 
Senetration flux, and average temperature af gases, 


respectively, on the air flow factor. These figures were 


Fe 


Setermined empirically for a specific surface of 0.63 ft /I1b. 
Curves are given for three values of the epecific fire load: 


2 
Pecos), and i2.3 ib/ft . 


Fiqure 3.2 shows clearly that the maximum penetration 


mux anad fully developed fire duration occurs at € = 1. The 
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specific fire loading appears to have little effect on the 
maximum flux, and the duration is usually quite short. It may 
also be concluded from Figure 3.3 that the hottest fires occur 
in ventilation controlled compartments. It is passible at 
this time to draw some conclusions about shipboard fires. 
Below decks fires which da not communicate freely with the 
atmosphere are ventilation controlled. These fires are more 
destructive as the duration is greater, the average 
temperature is greater, and the effective heat flux is greater 


than for an above decks fire. 


McCarthy demonstrated the usefulness of Harmathy’s 
technique for predicting tire severity of a typical berthing 
Sempartmemt. Figures 3.2 and 3.3 can be used for a typical 
berthina compartment, and in accordance with Table 3.1, a 


Z 
specific fire load of 3.0 Ib/ft was used. He assumed that 


the ventilation is secured, and that one water tight door is 


opem for access to firefighters. 


Z 
Xe = 15.0 Ft 
Ww 
h = 10 ft 
WwW 
U = 11384.2 \b/hr 
a 
2 
A = 150 Ft 
F 
oO = 0.75 


Hence, the fire is fuel surface controlled. Notice that 


as A increase the fire approaches ventilation controlled. 
F 


meem Figures 3.2, and 3.3: 


3 
T = 1.135810 sec = 19.2 mini 


Fe 
q = 6344 BTU/hr—-Ft 


s) 
ie! clo FE 


It should be appreciated from this discussion that fire 


is tndéed a complex process with several unquantifiable 
Variables. The term destructive potential was Introduced in 
reference 1. The destructive potential of a ‘fire may be 


considered a function of the maximum temperature obtained, the 
duration of the fire, the duration otf time at maximum 
temperature, and the quantity of energy transmitted from the 
fire. Most Fire Engineers use 2a time temperature curve as a 
definition of heat exposure. For the purpose of this stud», 
the use of an accepted time-temperature curve will be 


considered an adequate expression of the fire process. 
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3.2 Development of a Time Tempecatuce Curve 


Traditionally, fire tests have been performed based upon 
prior established standard time temperature curves. The most 
widely used fire test in this country is the ASTM-E119 which 
Was developed by the American National Standards Institute, 
(ref. 3). This method was established as a result of the 
necessity to standardize the structural performance of 


building construction. 


The basic procedure of any fire test if ta place the 
element to be tested into a gas furnace in which the standard 
time temperature curve is duplicated. The gas temperature is 
closely monitored through the Wuse of thermocouples. 
Additionally, deflection, axial thrust, and collapse behavior 
are measured. The E11? time temperature curve initially» rises 


e) 
Papidly reaching 1000 F at S minutes into the test. Beyond 


this point, the temperature is increased more slowly so that 


O 
at the end of one hour the furnace temperature is 1700 F. The 


a) 
maximum test temperature, 2300 F is reached at 3 hours and 


held constant thereafter. 


Another test method commonly used by Ship Structural 





Engineers i8 the International Convention for the Safety af 


Life at Sea SOLAS) ASO test. This test is used for testing 


non-structural Sulkheads from berthing and Gftice 
compartments. The specimen observed is not loaded but 
strained at the edges. Other non-structural bulkhead tests 


are SOLAS ASU, AIS, AO, BIS, BO, But the requirements for 


these tests are not as severe as for the AdO test, (ref. 1). 


For the purposes of this study, reterence 6 supports the 


oO 
use of a step change to 1700 F vice the gradually increasing 


temperatures of ASTM-E11? or the SOLAS AS6O curves. Figure 3.4 
gives the ASTM-E11? curve, the SOLAS ASO curve, and the 
assumed time temperature relation to be used if this 


analysis. 


| 
0 
On 

{ 





KPNUSCC 


7 


1 ee yp Ss 


mck 
[el 7— ee 
- aie ——— — 
ss  — 
a y ere 


[/ ere een asa A.C! 


—— — ASTM-E119-76 
A 
—- 4935 UM=z5 Curve 


£) 
(y' 


aie 





~]| 
Ww) 
~ 


*® an 
wae. 


o’ 
4) 
ae 


W 
~J 
; - 
a | ee | -—-~--}- 
ee 


a, 
fx yr 
‘ 


a5 
O 
Oj 
O 
S 


Zo =e) 


fo) 


Nae, min 
\ime/ \EMP ERATURE LURE FOR 


\/ARIOUS \IMPES CF Fizes 


- Ea SR 


Chapter 4 


HEAL TRANSFER AN4L YSIS 


The intention of this analysis is to determine’ the 
relation Between temperature and the distance along the length 
of the beam. A beam exposed to a fire Joad will be modelled 
as a fin with a temperature gradient. The temperature at the 


exposed face is T 3; the temperature at the end of the fin is 
b 


T , the ambient temperature. The geometry of the fin is 
a 


Sertined in Fiqure 4.1 


$.1-Conduction Heat ITransfec 


The following assumptions are used: 


i 
a) 


- T = Temperature at xX 


i 
=i 
if 
i 
- 


Temperature at x 


- The heat transfer is steady state. 


1 
0 
ee) 

{ 





Tin CoONFIGUR@SNON 


Figuem “4.1 


- K = Uniform conductivity of the material, 
QO 
Bev itin-sec-— Fi 
- h = Uniform surface coefficient of heat transfer, 


2o 


eer thra=tt - Fd. 


Additionally, this study will use the assumption that the 
rod is thin, which implies that the temperature in the thin 
direction, perpendicular ta x, is uniform. Therefore, the 
surface temperature equals the centerline temperature, (ref. 
7). APPlyIng an energy balance on an elemental section *x, 
yields a governing ditterential equation. In steady state, 


with no internal heat generation: 


Gece SAA Cao) 
Expanding the first term of the right hand side gives the 


following equation: 


dq 

_ = ~ as 
7” Gtx — & 7 (4, + 3) (4.2) 
Mx” Let 7 “Bx Ox (4.3) 
Gy 7 Ft Ax 7 oll Ea ie SEs Mie (4.4) 

Ox OX 
q. 7 = 98 pg Ot 1G BY 
De Aeneas) Ac ; 
ox Ox 

Where S = Surface area = AyAz. Conduction heat transfer across 


AA it given in the following 
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eeuation: 


Ze 
AL = = (9.6) 
a 
= a C4.7) 
or onl ak TIA 
2 ec hPAx(T - tT.) (4.9) 
Where 
P = Perimeter = 2¢(Ay + Az) 
T = Temperature of fin at x 


Substituting 4.35) and (4.8) into (4.1) vields: 


d oT (4.7) 
ae - T -T = 0 
zp (kSx) Ax hPAx ( A 
Finally, the governing differential equation for g@= (T - 
a 6} SS 
a 
a9 hP 
SEL ae eae s) = 0 , 
Sere S, K, and T are constant. 
a 
The differential equation far one dimensional steady 


state heat flow with one internal heat generation source is: 


Gia = Gani 


Note that the fin problem is similar to the single internal 


heat 
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generation source problem with a negative heat source: 


_ AP 
Wer 4k 


fee Se Lub op 


The methad af salutian will follow a method 


reterence 8. Assume the solution is of the 
equatian (4.13): 
9 = een + ee 
Where: 
B = [hp/ks}]?/? 
The boundary conditions for a fin of finite length, 


defined by Figure 4.1 are: 


@ x 


Ht 
cp 


O; 0) = 


@ x L; oh = Whe os Xa) 


Substituting (4.15) and (4.16) into ¢€4.13) yields: 


- G4? - 


(oars. 


provided in 


form given in 


(494.13) 


€4.14) 


L, as 


(4.15) 


€4.146) 


C4.1 9) 


(4.18) 





eying tor C 3: 
1 


Si = b CO.17) 
Bea eo 
S —C 


mee tollows from (49.17). Finally the solution is: 
Z 


g  gbih-x) , ,-B(L-x) (4.20) 


G oBb -BL 


At this point, it is mecessary to solve for B. For the I 


beam in Figure 2.3, the perimeter ta surface area ratio is: 


Pee (Oe Za t) 

Cee ae €4d.21> 
Oe ea2e ts Ee 2k a 
W £ow sig 


Figure 4.2 gives the relation between thermal conductivity and 
temperature for plain carbon steel, (ref. %)2. Within the range 


QO a) 
af interest, 70 F <=-T <= 2000 F, the curve is considered 


approximately constant. For the purposes of this analysis, 
O 
the thermal conductivity at 1000 F will be assumed. 
k = 0.0005 BTU/in sec °F CA ae: 
The condition of nmatural convection in a gas iS assumed 


mer tmis analysis, (ret 7): 
h = 2.0 BTU/hr ft- OF (qe) 


Therefore, 


B= (P/10.8s)*/7 (4.249) 


where, P/S is given in equation (4.21) 
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The validity of thin fin theory must now be checked. 


According to reference 7, the thin fin theory is valid if: 


ee 
aE 3 (4,25) 


lv 


Substituting equations (4.22) and (4.23) bate? €4.25).) the 


maximum value tor b is: 


< 
bt= 306 £6 
(4.26) 


This 18 a reasonable criteria for this study. Therefore thin 


fin theory is valid. 


mee kesulis 


It is now possible to draw some conclusions about the 
heat transfer through the beam. A Microsoft BASIC computer 
Program was written in order to examine several geometric 
configurations af the beam. The source code for the program 
is im Appendix A. The results are given aqraphically in Fiqure 


4.3. The dimensions of the beams examined are in Fiqure 4.4. 


The following conclusion may be drawn from Figure 4.3. As 
P/S, and L increase, the effect of the fire becomes more 


localized. With a typical shipboard I beam, the ratio of 
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perimeter to surtace area is such that the fire load may be 
treated as a point load. That is toa say that the heat 
transfer through an I beam remains highly lacalized in the 
mecinity of the fire. Al] effects of the fire virtually 
disappear at a distance which is a quarter of the length of 


the beam away from the fire. 
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Chapter 5S 


ANALYSIS OF THE Fell ure MECHANISM 


AS was stated in the preceding chapter, an I beam exposed 
to a typical shipboard fire will experience the same failure 


mechanisms as a beam under a point load. Recently there has 


been increasing interest in this type of failure due to 
mechanical loads. It is the intention of this study to 
approximate the shipboard fire load 4 an additional 
mechanical load, and to use two of the recently developed 


techniques for analyzing the failure mechanism. Additionally, 
this chapter will include summaries of literature which 


focuses upon bending collapse of thin walled structures. 


Kecman (ref. 10) studied the bending collapse behavior of 
rectangular and square section tubes both theoretically and 
experimentally. He used a limit analysis technique and a set 
of formulae to derive the relation between the hinge moment 
and the associated angle of rotation. Very good agreement was 
found between his theoretical predictions and the experimental 
results. His theoretical predictions were in the form of a 
mumerical solution. However, Hayduk and Wierzbicki (ref. 11), 


gave a closed form analytical solution. In reference 11, 
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lower and upper bound solutions for the mean crushing strength 
of cruciforms were obtained by considering bending and 


extensional modes of deformation. 


oai_Kecman”s Method 


Kecman developed a numerical scheme toa analyze the 
bending collapse of rectangular and square section tubes, 
(ref. 10). The primary failure mode is Jlacal buckling and 
hinge development. He considered only uniaxial bending 


collapse. 


31.1 Maximum Bending Strength 


Kecman referred ta a study which examined the bending 
strength of a box shaped section with the webBs in tension and 
the top flange in compression. The effect of the webs was 
found to be negligible. Therefore, the critical stress of = an 
I beam may be assumed to be the same as the critical stress of 


a box shape beam. 


The critical stress, go , of the section was derived ina 
Gr 
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form: 
2 


2 
T EF oe 
GC = (5.23 + 0.165) (=) 


a = flange width 


6b = web height 


t = web, and flange thickness 


If the critical stress is less than the yield stress, gg of 


a 
the material, then the flange buckles. The effective width 
becomes: 
Sor 
gO ee * Ons) (5.2) 
= 
oO = stress at the edge of the flange 


2 
If the critical stress is greater than or equal to the yvield 


stress, then a =a. For the purposes of this study, assume 
f 


that collapsing occurs when oa == 10458 tt tal lows that tie 
2 a 


maximum moment, M , may be calculated using (3.1) and (35.2). 


m 
ic o iS /acamee 
cr Cc 
a 
22+ b+t+a (3— + 2) 2 
m Oo 3(a + b) 
i+  C I= 3 Oo 
cr fa) 
a 2 (Capel whe, 
te Roe ee) tO (bee 2c) 


where M = the fully plastic yield moment of the section. 5s 


im < oa < 30, a moment, M “, is defined: 
o cr o a 


MeeG cote =) (5.30) 
O Oo 3 


With the resulting maximum moment: 


0 =" 6 
M=MsM! e ne O 3.30) 
m oO * (M. o ara 
O 


3.1.2 Actual Bending Collapse Mechanism 


A detailed analysis of the hinge deformation pattern 
offers insight into the failure mechanism. The deformed 


Becterh is given in Figure S.1. 


The following characteristics of the deformation 


mechanism were given in reference 1: 


Cad Most of the plastic deformation is concentrated along 
Meo taticmnarye yield lines: JF, HI, EB, GS, FC, HE, CR, 


Pe, GK, and EL. 


Cb) Travelling yield lines of rolling deformation were 
also observed along GA, AE, KA, and LA while the distance 


BA changes during hinge rotation. 


€c) The corner angle along KL remains almost constant. 
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Cd) The deformation of the tension flange includes the 
travelling yield lines NP, and MQ, approximately 
corresponding to the intersection of the undeformed 
flange planes and the current cylindrical surface TU. See 


miugures Oeil, and 3.2. 


Ce) In-plane deformation of the walls was also observed, 
particularly along the rolling lines GA, AE, BA, HC, CF, 
IR, and JR. This deformation is essential to maintain 


hinge Kinematics. 


Cf) The rolling radius, r, varies slightly along GA, 
measured By cutting at right angles to GA and its 
estimated value was approximately proportional to the 


hinge length. 


Ro =~ (obe0S-0 037 4KL 


The radius reduces during collapse, so that rolling 
becomes more and more difficult, until the mechanism jams 


and starts a new phase of collapse. 
(gq) The rolling radius along KA, and LA varies from 


infinity at K and L to a value close tor in (Cf) above. 


Four distinct phases in the hinge development were 


identified: 


C1) Initial phase with a bulge in the web, A , in Figure 
1 
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Dae 


(2) Second phase with rolling as described by the 


cross-section shape A in Figure 5.3. 
B 


(3) Third phase indicated by a jamming of the rolling and 


described By the cross section BAD in FigureS.3. 
=) 2S) 


(4) Final phase initiated by the contact of the two 
buckled halves of the compression flange. There is no 


clear transition from one phase to another. 


3.1.3 Theoretical Bending Collapse Mechanism 


Since structures most commonly experience hinge collapse 
described by the second phase which includes rotation angles 
in the range, S-10 <= @ €= 25-35. Kecman based his theoretical 


analysis on the following assumptions: 


- The web and the flanges of the section deform along the 


concentrated yield lines only; 


- The web and flanges are incompressible and 


inextensible; 


- The structural continuity is maintained in two 


characteristic sections; and 
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- Some in plane deformation is essential for hinge 


Kinematics. 


The model is given in Figure 5.3. 
tre rolling radius is negligible in comparison with the 
other hinge dimensions. Points A, B, C, and R are in a 
straight line. The coordinates of point B are: 
. . Wye 
B(x,y,Z2) = (h, beosp - [bsinp(2h- bsinp)] SOG) (5.49) 


From continuity of the middle section, the web height, b, is 


determined: 


7 (ye 4 2,1/2 z 
lo, ee Ge Yn an Ceey, 
Pane sy = ¥ , and equations (5.49, and (3.5): 
= E 
ye = 
Zz. = psinco™ WSino™®’ (being (@n - bsinl¥’’ “*cosp eer 
An analysis of the rotation of the fiber in the top 


flange will lead to the equation for the hinge length, 2h. For 
a fiber which is initially parallel to the beam and passes 
through point A, IJlongitudinal continuity is implied. The 
distance from this fiber and the edge KL is AD=5 - ae 
Initially, the coordinates of the fiber and the line GK are 
Peo ateee emitter fotation, the coordinates are [C65 - 


Potion 2 )COsp, UI. Point A has the coordinates (Ch, » , 
- 4 A 
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mei. the distance from paint A and the rotated point :!s h if 
A 


the web and flanges are inextensible. The equation for h is 


therefore: 


Z : 2 2 2 
h = (h- (b-z,) sinp] - [y,- (b-z,) cosp] + Z, Coat 2 


Substitution of (35.43, and (5.6) into €5.7)3 and meeting the 
condition required By (5.5), itt may be said for any ratio a/b, 


and any rotation angle: 
2h =a or 2h=b Co. 


The value which yields the lowest energy absorption will be 


used. Thus, 


NWA 


IA 


2h = b PENS 


The second phase of detormation is complete when the 


rotation angle equals the jamming angle, a . Jamming starts 


when the two Buckled halves of the compression tlange contact 


each other. From Figure 3.3 the jamming angle is: 


a . Shes 
oe = 2arcsinp (= *) : Oe 


3.1.4 Nominal Energy Absorbed During Hinge Rotation 
The assumptions listed in the previous subsection imply 
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that all the energy !ts plastic and is absorbed along the yield 
lines of the theoretical collapse mechanisms. Relative 
rotation along the flange edges JF, and IH results from hinge 


meracion, see Figure 3.1. From geometry, of Figure 5.3: 


leat 
- 0 - arcsin(l - peine) (5.11) 


2 
i 
NT | 


The energy absorbed along these edges becomes: 


TT : b 
= = ee = - Coal e. 
W, a ees am als 0 arcsin(l psino) | 2 
m = Average fully plastic moment of a unit width of 
p 
the flange 


Similarly, itt may Be said: 


b ae ee, 
Wi, = V = = = ol 
5 ae eae 2arcsin(l ino) ] 


The energy absorbed by the yield ]ine AB can be found 


from €5.119, where the length of AB = z is Known from (35.4), 


a 


and the relative rotation is €7 —- 26). Hence: 


2 
/ SSE) OL On 


ee . 
W.=W = m (bsin o=inedunje-: (ada (Tava h 


3 AB 


similarly, along 66&, the length is h, and the relative 


rotation ts 2: 


W = W + W lie elxel Oi 


Figure 3.4 shows the relative rotation, "7, along yield 
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line GK. The construction of this figure is detailed as 


follows: the plane w is normal to GK and contains point A. 
xy 
Point A” is the projection of A onto the xy plane. A is 
o 
obtained by rotating wabout its intersection w with the xy 
XY 


plane. Coordinates of A" are: 


htanp + beosp - [bsino (2h-bsino) 1/7 
= f 
Yau aan oe (3.16) 
- ie ee 
ante Y,utanp 


The energy absorbed along both the GK, and EL lines are then: 


Z 
W. = W + W = eo ee ee, €9.18) 


[ (n= 4) a (Yaur¥_) 11” 


The rolling radius along the yield Jines GA, and AE 


remains relatively constant, so: 


Wo = W.. + W._ = 2m eve Reig 
Jo ae A 


During hinge rotation, the rolling radius varies as 


described by: 


r= (0.07 - Sa)h CO eC 
Equation Co. 20> is based on empirical data. Kecman 
approximated the variation of the rolling radius along KA with 
cane oermeolic approximations: 

ea 7 TP (5.21) 
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Where | is the distance trom K along KA. The linear variation 
K 


of the ralled length Between K and A is given by: 


i 
Q Dye la Conroe) 
r KA A 
The energy alang KA is: 
2m zZ_KA 
wy pA Gone eo) 
KA £18 
or 
2 ea 2 omy 
W. = a eS - 
7 Ramen ay th NGS + 2. ] (5.249) 


iieerelaciv¥iesrotation alemg yield lines PN, and QM ise, 


whereas along KL, and MN the angle is given by: 


s “1 
ee eeeareccan(—) 
A 
Srnaally, the energy absorbed by these yield lines is 


approximated By: 
Zz 


= emilee + 2harctan (—)] (9.25) 
V 


W. = W +W_ +W +W 
8 PN OM KL MN 
-A 


At this point, the total energy absorbed by the hinge may 
Be expressed as: 


: F — 
LC) =) W, (a) Gane) 


i=l 


This result is valid for a rotational angle less than the 


Jamming anqle. Kecman determines the moment-angle 
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relationship numerically from the following equation: 


M(q) = What fa) = Wa) Cale ie 
AG 
Where \@ represents a small But finite increment of the hinge 
mele, GW. 
This formulation assumes that the material is ductile 


enough not to fracture during hinge ratation or subsequent 
bending collapse. This iS considered a Safe assumption for 
common structural steels. The average plastic moment is given 


by the tollowing equatian: 





mee nC 
MD ole (3.29) 
Z | 
The yield stress, and the maximum nominal tlow stress as 


functions of temperature were determined in Chapter 5S. 
fet. > The Final Form of the Hinge Moment 


The analysis in the preceding section is not valid for 
the first phase of detormatioan. This phase will be 
approximated By assuming that for a moment aqreater than the 


aut Dlastic moment, M , is constant and equal toM . This 
O a) 


would be a conservative assumption. Far hinge angles qreater 


than the jyamming angle, Kecman offers the tollowing formula: 


M(a) = M(a.7) + 1,4[M - M(a,) J (a- a5) 
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This formulation has been programmed in a Microsoaoft-BASIC 
program. The source code is given in Appendix B. Figure 5.5 


shows the results of the program. 


5..2-Hayvduk.s_andWierezbicki~s Method 


The primary difference Between Kecman’s method and this 
method is that HaydukK and Wierzbicki considered mot only 
inextensional deformatioan modes, but also extensional 
deformation patterns. The effect of these extensianal mades 


will be fully developed in this section. 


Hayduk’s and Wierzbicki’s formulation results in a closed 
form analytical solution which predicts the mean crushing 
ferce. Reference 11 specifically addresses the problem of 
analyzing the crushing process of thin walled, plate formed, 
open structures with particular emphasis on L shapes and 
crucitorms. The I beam can be approximated by two T shapes or 
four L shapes. Lower and upper tboaund egolutions far the mean 
crushing strength of cruciforms was obtained bY considering 
various modes at deformatioan. According ta Hayduk and 
Wierzbicki, extensional deformation accounts for at least ane 


third of the dissipated energy. 
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The end products of this section will include the 
application of Hayduk’s and Wierzbicki’s method to an I. beam, 


and reviews of both references 11 and iz. 


2.2.1 Collapse Modes 


Hayduk and Wierzbicki considered out or plane 
deformations in which the deflections exceed 10-100 times the 
section thickness and are of the order of magnitude af the web 
height and the flange width. The process of collapse starts 
as the compression flange loses stability, either elastically 


Seep iasticaliy at a critical load, P . Following the buckling 
E 


deformation, the flange can maintain added loads until a 


maximum load, P » 163 reached. Eventually, as the resistance 
max 


decreases, the structure collapses completely. The structures 
being examined here often regain strength after the initial 
loss of resistance. Subsequent load peaks reveal the 


development of additional folds. 
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S.2.1.-1 Isometric Transformation Modes (Inextensioanal!l 


Deformations) 


These moces are of the type used By Kecman to analyze the 


collapse mechanism. In this type of deformation, the energy 


is concentrated over hinge Jines, while the rest of the 
Structure remains rigid. An assumption of inextensicanal 
deformation is that material continuity is maintained. Ouring 
the transformation of two surfaces, the lengths of 
corresponding arc elements is preserved. Hayduk and 


Wierzbicki show the progression of various deformation paths. 


These are illustrated in Figure 3.6. 


Timp i aine os.gbca)yeana <b), the final state is formed 
inextensiBly. This can only be done if the boundaries have 
Sea ricient treedom to deform. The element in Figure S.6 (Cc) 
is initially folded along the stationary hinge line AB; 


finally the element is rotated extensibly, ‘ref. 11). 


g9.2.-!.2 Hinge Collapse Modes 


There are two types of hinge collapse modes. The first 
is a propagating hinge line which occurs when the hinge line 
is free to move with respect to material points, and the 
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material of one flange is effectively transmitted to another 
flange. Extensional deformation occurs in the vicinity of the 
hinge line. This mode is a characteristic of thin walled 
structures consisting ot two intersecting plates such as the I 


beam being considered in this study. 


The second hinge collapse mode is the formation of a 


Stationary hinge line. This mode is characterized by a hinge 


line fixed in the material. Side tlanges which are initially 
rectangular are transformed into trapazoidal elements. 
Therefore, extension occurs over an entire flange. This made 
will not be considered in this anaiysis, (ref. 11). 


9.2.1.3 Assembled Collapse Modes 


Hayduk and Wierzbicki demoanstrated the methed of assembling 


'D 


the folding patterns shown in Figure 3.6. The patterns ar 
folded into more complicated symmetric, asymmetric and mixed 


collapse modes which are shown in Figure 3.?. 
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5.2.2 Analysis 


For the purposes of this study, the analysis of aT 
section will be used; this analysis is considered to be 
equivalent to an I beam analysis. In reference 12, the short 
T cross section under consideration Was subjected to 
compression with clamped boundary coanditians. A model is 


presented in Figure 5.8. 


For the purpose of estimating the mean crushing force, 


P , the T section is expanded into a plate and an angle as 
m 


shown in Figure 3.9%. 


The plate will collapse through the bending zane limited 
by the two horizontal moving lines. The angle will follow the 
folding mechanisms described by Figure 3.10¢6b). The collapse 
modes in Figure 93.10¢6) is cansistent geometrically and 


Kinematically. The regions are described as follows: 
(I> Four plane trapazoaidal elements as rigid bodies. 


CII> Two sections af cylindrical surfaces at which 


continuous Bending takes place without any extension. 


CIII> Two sections of conical surfaces in which material 
is Dent and rebent again as the material moves fram ane 


flange toa another. 
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Soy enesee ci eaneot taoroeidal surface which produces 


extension in a circumfrential direction and continuous! y 
changes principal curvature in the other direction, ‘ref. 


11) 


The following assumptions will be used in this analysis: 
- Neglect the short side of the angle element. 


- The length of the two horizontal hinge lines is 
approximately equal to the length of the large side of 


the angle element. 


- Two inclined hinge ]ines are of a length equal to a 


half wave length, H. 


- The length of the toroidal surface in the vertical 


direction is neglected. 


The energy dissipated of an angle element it given by 


summation af three terms, Cref 11): 


1. The energy of a continucus deformation field in the 


section af a toroidal shells: 


EF = 9. 28M Hr/t (5.230) 


2. The energy dissipated by the discontinuous velocity 
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Poth the Mortizental hinge line oft length, 1: 


3. The work done by an inclined hinge line through the 


Geformation process: 


Z 
E. a 2.22M H Sa (5,32) 
where: 
r = The small radius of the toroidal surface. 
2 
M = oft J £4, the fully plastic rield moment per unit 
O Oo WwW 
Width 


H = Mode halt-wave length 


TF = Oa7o 
oO Ou 

a = Ultimate strength 
Ou 

t = Web thickness 

Ww 


For the models given here: 


ae ar Be web is 4(E .) evange © 3 CaeRee 


eifter Scuestitution of the dissipated energies into (3.33): 
Z 


1 oe 
P/M, = 4.642 + 2c + a() 1 + 2.228 
- H iz xr 

WwW 
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OBtimizing equation (3.34) By using the 
derivatives: 

3m _ 0 

or 

atm _ 

oH 
Beating ©€S.35a) and (5.355D) yields: 

eel, 3 
——— ame Te 7 (c + ayes )] 
4.64 be 
2 


ie 


Saree 


Ultimately, the 


Gissipation appear equal 


mean crushing force: 


ee = 3[(4 


mauation (5.38) may be si 


Pp = 
ae ae 


nneres 


Note also that equations 


ae Se 
eammcay oe I) ”) 
WwW 


three 


a 2,1/3 


basic mechanisms of 


ly weighted in the final 


iz Z 
3 
64) (2.22) 1(1 + k(S) ) (Ly 1°7 
WwW W 
molified as follows: 
Cathe 3 
SIAT ADA, |e 
W 
ome 
c 
A, = 4.64 
t. 2 
A, = T[{l + k (=) ] 
W 
ae = Deo 
Ges ieessse>, and (5.37) 


following partial 


(54.305) 


Ca. Say 


energy 


equation for 


(5.90) 


C3.41) 


may be 





eemoiitied into the following forms: 


_— a cS H = 
eee 7 le if AO ? ae ia? 
Z Pe DS (3.45) 
1 [(A,A,/A, y(ct )] 


az L 2 1/3 @aeicy 
H = 
[(A,/A, A) (c De 

The application of these equations to the case of an I 
beam requires that the geometric parameters be madified. An I 


beam is two T sections, hence the dimensions are daubleds; but, 


b = 2c, which yields the fallowing relations. 
" 3 
pe7 = 3 {AJ A,A, (b/2t )] (5.487) 
he = AAS (3,400 


It 13S now possible to apply this relation determined for 
the crushing case to the bending case. The geometry of this 


epuching mode is given tm Figures 3.11 and 3.12. 


Recall the propagating hinge line mode as described in 
paragraph 3.2.1.2. This mode describes the global geometry of 
the crushing process and is determined purely by the width of 
the two side panels, d, which is prescribed, the wavelength, 
H, which is to be determined, and the crushing rotation angle, 
QA, which 1S &@ parameter of the process. The relation Between 


the axial displacement, 5, and the angle, @, is: 


6 = 2H(1 - cos6) (S249) 
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D 


Crusuing Case 


Figure 5.4\ 





Extensions. Co.rapseE Mone - 


Derin TION OCF GEOMETRIC LA RAMETERS 


Wel ere 





Reference 12 determined an expression for the instantaneacus 


crushing force, P, normalized with respect to P 
m 


al 2 0.9cosR (9.90) 
Ee 3 2cos8 Ye fanoil + sin’B) 177! 





Fiqure 3.13 gives the geometry of the bending case. From 


mroqure 3.13: 


veges 

b Cae ales 
mus, it may be said thats 
p/p, = £(6) = £(2btan5) (5.52) 
Potro a section width, bs 

M = = . (5.53) 

P 
M =—b (5.54) 

m 2 
Where M is the mean moment due to bending rotatian, a, as 


mM 


defined in Figure 5.149. After equating (5.53), and (5.54): 


C9.09) 


sa 
it 
a 


Figure 5.15 gives P/P as a function of 6/2H based = on 
m 


Seon eG. o.oo), €rer., IZ). This oraph is an 


equivalent description at M/M as a function Gr 
m 


CO/HITtanta/se2d]. 
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SENDING CASE 


mouree 5.13 


agi Wied £ 


Verinition of Angle of Rotarion 


Ficguee 5.14 





te 


NN 


ae 
LS io ‘a KS) 1.0 


Pe | 4 
fs = 
Cd » 


Mom ENT Disriac EMENT RELATION 


Fiouree 5.15 





For this purposes of this study, an approximate relation 


for M/M anda will be determined. Assume equation (5.54) is 
m 


a valid function for M/M : 


m 
M/M_ = D/ (6/2H)® (9.98) 
where: 
Db, om are Bere ence 
Assume the new curve goes through paints P , and io in Figure 
1 
alo . 
P,(6/2H, M/M_) = (0.30, 1.0) (3.97) 
P,(6/2H, M/M_) = (1.0, 0.24) (9.98) 
peiving for D, n: 
D = 0.24 Case. 
a 5 Cano) 
Assuming n = 1.0, the following equation results: 
M/M_ = = Coeol 2 


Figure 3.15 shows how the approximate salution compares to the 
Sete Chelan. oust i tuting equation «35.513 imto equation 


nol onicrol lowing results: 
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IS 


3H 
ComeaRison of AperoximMaTe 
Ano Exact Sow Tons 


TiGURE S09 


0.24 


te = (b/H) tan (G/2) 


Boensidering equation (35.54): 


wa Te ooze 
2 (b/H) tan (a/2) 


Now substitute equation (35.47) into (35.63) 


ON essay (es 


M= [3M [A,A.A_(b/2t ae 
O Ww bteam(@/ 2) 


tere 


But M is the fully plastic yield moment for the flange 





The desired final result of this analysis 


between M and M ”. Rewriting equation (35.65): 
a) 


M'=o f zdaA 
oO od 


Reeail eguations (2.2) and (2.4): 


eee) = =20-5F + 40,000.00 


at zdA = abt + 


Diittivnwecoiipstt tite equation (2.4), (3.46), 


Soto) hice Co. o4) : 
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Cane) 


(3.63) 


(3.649) 


Se a) 


relation 


(3.66) 


Sin era and 


<10o-" 


Whi 


EO 


ray 1s 





: “a / ; 


te, 


a eae ee: se 
NN 


| 


Fae 


PS OO kf | |S SG 
K Radians 
ecu, OF 4 Ay~DUK. AND 


WMieRZBICKIS ANALYSIS 
Vee 5.17 


Cl) ae 


rr 


O 


30, t* x 0.24[(AS/A,A,)b-2t]~/- 
wi IS.) Sie a aeons Ne 
ee 2 [A] A,A, (b/2t) J (>) | btan (Q/2) 


40 (abt + b°t/4) 


The simplitied final result of this relation iss 


0O.367t(b + 2a) (5.68) 


nN i. ftan(@7 2th 4+ 4a) ) 


Figure 5.17 shows araphically the resulting moment-angle 
rotation relationship as determined by Hayduk and Wierzbicki. 
Appendix C contains the source code for a Microsoft BASIC 
program which determines the moments resulting from equation 


oa 2) 


Boe Jelemnoat Son OtlPe sults of the dwo Methods 


The primary difference between the two methods oat 
analyses is that Kecman considered only inextensional 
deformation modes, whereas HaydukK and Wierzbicki considered 
extensional deformations as we] | as inextensional 
deformatian. Since the second methad considered more energy 
GSiescipation —than the first, the allowable moments due to 


detormation are greater. However, for angles greater than ane 
tenth ot the jamming angle, Kecman’s method yields greater 


moments. The errors induced by his mumerical integration 
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technique are compounded as the angles increase. Figure 3.18 
shows graphically how the results of each analysis differs. 


The dimensions of the section analyzed are given below: 


a= 12 inches 
b= {2 inches 
t = t = .6 inches 
f W 


The maximum moment given in Figure 3.18 ts the fully plastic 


yield moment of the section, hereafter called M , which is 
rw) 


only a function of geometry and temperature of exposure of the 


beam. Hence, M i185 not a function of the analysis method. 
fm) 
Clearly, as the angle of rotation increases the moment 
decreases; and as the temperature increases, the moment 


decreases. Appendices B and C also contain results which show 
the effects of changing the geometry of the section. For both 
analyses, the moments are decreased as the section thickness 
decreases; the web height decreases; or the flange width 


decreases. 


These methods of analyses are considered accurate 
predictions of the tailure mechanisms of a beam exposed toa 
i Meee Recommendations far using the information provided by 


these analyses will be given in the following chapter. 


- ¥j{ - 





We 







ta Keeman's Me-tiep 


Tay Come 


i \ i 
eee MArDuUK e WierRz B\CK\S 
~ Mertuoo 


O 
N 
(d 
pe oe 
a 
() 
J 
0 ae 
£) 
- 


o, RADIANS 
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Chapter 6 


CONCLUSIONS _ 6ND_ RECOMMENDS TIONS 


A summary of each chapter is presented here with 
concluding remarks. Finally, recammendations are affered for 


using this analysis for future studies. 


4.1_Summary 


As etated in the introduction, the ultimate objective of 
this study was to show how much the strength of a shipboard 
structural member is diminished when exposed to fire. To this 
end, Chapters 2, 3, and 4 provided input into the analyses of 
the failure mechanisms studied in Chapter 3S. The results of 


the preceding chapters will be summarized here. 


Chapter 2 developed the model! of the structure. Simple 
observations revealed that the ship is a redundant structure 
composed of beams, plates, and columns. The motivation 


towards choosing an I beam for the analysis was developed. 
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The material selected for the model was ASTM-ASS steel as its 
physical and mechanical properties are similar ta the 
properties of ABS Grade A steel, a common shipbuilding steel. 
The equations for the temperature dependent mechanical 
properties were the output of this chapter to be used in the 


analysis. 


Chapter = provided background oils characterizing 
shipboard fires. Three parameters which are assumed to be 
Peeumanvesoregiections of the contributing factors of a fully 
Geveloped fire were defined. These factors ares: the duration 
of the fully developed fire; the overall penetration flux; and 
the average fire temperature. Since fire is indeed a complex 
process with several unquantifiable variables, a time 
temperature curve was used as a definition of the fire 
process. For the purposes of this study, a step function toa 


s) 
1700 F was used as the time temperature relation. 


Chapter 4 was the heat transfer analysis which determined 
the relationship Between the elevated temperature exposure and 
the distance along the length of the beam. The results of 
thie analysis led to the conclusion that a beam exposed to 
fire will experience the same failure mechanisms as a beam 


under a paint load. 


Finally, the primary thrust of this study was to analyze 
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the failure mechanisms due to fire exposure. Two methods were 
fully developed in the preceding chapter and the results were 


slmmarized graphically in Figure 3.18. 


46.2 Recommendations 


The ultimate objective of any structural desiqn pracess 


is to develop a structure which is strong enough to withstand 


all expected service laads and a high percentage of unexpected 
loads. The information provided By the analysis in Chapter 5 
is considered usetul input far the desiqn of shipboard 


structures. It ts possible through the use of Figure 3.18 to 


determine how the stiftness of a structural member is 
decreased when exposed to elevated temperatures. The 
diminished stiffnesses may be used as input tao ae *#Ftinite 


element program which simulates a ship structure. The effect 
of a shipboard fire upon the entire structure may then be 


determined. 


The examination of the effect that one weakened member 
has on the other members of a redundant structure is 
considered fertile ground for future studies. Reference 16 
may provide insight into the behavior of this phenomenon. 


ZuBbelewicz and Mroz numerically simulated rock burst processes 
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treated as problems otf dynamic instability. According toa 
reference 16, local plastic failure occurs at one location 
(such as in a member of a ship’s structured and the elastic 
energy flows to the hinge and is the driving force into the 
damage. Energy then flows into the neighboring hinges, which 
is to say that one unstable member may produce a chain 


reaction and other members may subsequently fail. 


As a final word, it is highly recommended to utilize the 
information provided in this study ase an aid for examining the 
Parect Of fires on the shie”’s structure. Ultimately a ship 


may be desiqned ta withstand shipboard fires. 
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Appendix A 


HESatTReaNSEER eNelLysis_PRDOGREM 


REM PROGRAM TO CALCULATE TEMPERATURE GRADIENT IN &A 


DIM TEMP¢ 20) ,Y(20) ,2¢20) ,TOP< 20) 


Piet rMircR THE GEAM LENGTH IN FEET" 3A4L 


INPUT "IS BEAM SOLID ¢€$S) OR I (I)" ;NS 


IF N$="I" THEN GOTO 90 
LPRINT " SOLID BEAM" 

INPUT "ENTER WIDTH OF BEAM" 3W 
INPUT "ENTER DEPTH OF SEAM";D 
LPRINT " WIDTH, DEPTH="3 W;D 
B=SGR((12/5.4) #¢C1/DI +C1/W))) 
GOTO 155 

LPRINT " I BEAM" 

INPUT "ENTER WIDTH OF FLANGE" ;W 


ewe enter DEPTH GF WEB" ;D 


LPRINT “ WIDTH OF FLANGE, DEPTH OF WEB="; 


INF UR “ENTER FHICKNESS OF FLANGE" ;T! 


Pie ENTER TRICKNESS OF WEB" ;T2 


wee imlenNe ss OF PEANGE, THICKNESS OF WEB="5T1;T2 


P=2*(D+2*W+2*T1-T2) 
SS=Z2*¥CTIL¥WI+0CD*T 2) 
B=SQRC012/35.4) CO P/S) >) 

Ci eeeNorH OF BEAR" sx; " 


LPRINT " B=";8 : LPRINT " * 
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FEET" 


lJ 3D 





160 


La 


BL=B*XL 

SGT—EAP* SL)—-EXF< =—BL> 

X=XL/20 

EeEinn oo) eh INO XA T-Ta/io-Ta" 
FOR K“=1 TO 20 

YOCKADEX¥KX 

£0 KAI=BECXL-YCKAD 

TORI Avex S62 CKA3)-EXPC—2 (KZ) 
TEMPO KAI=TOPRPCK4D /BOT 

Perein K, y¥ © KAD, TEMPC KAD 
NEXT Kx 


END 


mee (lO 10 oe 





Appendix B 


KECMAN” S_aNelysis_ PROGESM 


10 DIM W¢20) 
20 INPUT “TEMPERATURE OF EXPOSURE, DEGREES F";T 
30 INPUT “FLANGE WIDTH, INCHES";A 

40 INPUT "WEB HEIGHT, INCHES";B8 

41 INPUT “THICKNESS, INCHES";TH 

42 LPRINT " ° 

SO LPRINT " " 

Sit LPRINT " TEMPERATURE=" ;T 

52 LPRINT " FLANGE WIDTH=";4 

53 LPRINT " WEB HEIGHT=";B 

54 LPRINT " THICKNESS=";TH 

60 S$O=-20.5*T+40000.0 

65 IF T>=1951 THEN $0=4.5 

67 MO= SO*B*TH*(A+B/4) 

68 LPRINT " FULLY PLASTIC MOMENT, INCH-POUNDS" ;MO 
70 J=0 

90 U=3.14159265 

100 IF & <= B THEN H=4/2 

110 IF A > B THEN H=B8/2 

120 QL=(H-.5*TH)/B 

130 D=2¥(ATN(QI/SOR(-Q1*Q1+1))) 

131 LPRINT “ JAMMING ANGLE, RADIANS";D 


132 LPRINT " ANGLE; MOMENT/MO; MOMENT" 
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X=SU4CTH* 2074 

pee K=0710 TO 112#b/710 STEP D/10 
Pi = K/2-K*3/48+K*5/38490 

P2 = {-K*2/8+K* 4/384 

eer Sar! 

P4 = (28H - PS) 

PS = BeP2 

P6é = H*¥P1 

pre = PI/SP2 


Wee— ero CoOR (PSP 4) > 

R=¢( .07-K/70)*H 
2A=(PS¥P1)-P6+¢(SQR(PS*¥P4) 0 *P2 
Qi=1-(B/H) *Pi 

QZ=ATNC QI /SQR¢-GQ1#Q1+19) 
WI=2tA*86U/2-K/2-G2) 

W2=A* €(U-2*Q2) 

WES 2 Sled= | See SIR) g86 Cr ts) dere 4) 0 lian 1 ae. 
W4=H*U 

eGo eet a. oon. PSFP4 I) I/ CL +tP7*2) 
LONG ES = 7 

GTC GC. HAC) 720+ ¢ CYC-YB)*2) 3) 
Q3S=ATNC QI) 

WS=2*B#Q3 

Wie, Sxteee hs SOR. CH 2) +0 ¥B* 29+ 2ZA*2))) 
W6=2*CH/RI*ZA 

G1i=ZA/YB 


Q3=ATNCQ1 > 


ime (0) a 


410 


420 


WS=22CAECK/2) + 2*H*Q3) 


J=J+1 


We JIHW1 tW2tW3tWl4+WSt+WstW7tws 


NEXT K 


eek l/= Ged =1 


Y=C(TH/ CBC S9#AtB DD RECWOCT AFIT AWOIT ADO OD/10) 


M1L=Y MO 


EPR iN: 


NEXT 14 


LPRaNT ™ 


bERINT ” 


LPRONT ° 


le 


END 


"sT%*D/103:YIM! 





10 


30 


40 


30 


61 


62 


70 


80 


70 


100 


) aif 


INPUT 


INPUT 


INPUT 


INPUT 


Pi LING 


LPRINT 


LPRINT 


LPRINT 


LPRINT 


Appendix C 


HSYDUK_ SsNDO_WLERZBICKI“S_SsNélLysSis_FPROGESM 


feVeerenURE OF EXPOSURE, DEGREES F“;T 
"FLANGE WIDTH, INCHES" ;A 
PeEESRETGH), INCHES" ;6 


Mim onNeESS, NGCHES* ;TH 


Se heMPERATURE ";T 
" FLANGE WIDTH" ;A 


" WEB HEIGHT";B 


EPR Nie. UHlCKNESS* ;TH 


T1I=(R*B*#TH+B*¥B*£TH/4) 


T2=-20 .5*T+490000 


eee) es) THEN T2=4.5 


MOQ=T2*T1 


Peete y Pest ic MOMENT, INCH-POUNDS" ;MO0 


T3=0 .36*43.1491593*TH*(Bt+2#AD 


T4=1/ 


(BHC Bt+4#A) ) 


LPRINT" ANGLE; MOMENT/MO; MOMENT " 


IF AC=B THEN J=A/2 


IF A>B THEN J=B/2 


Q1l=(J-.S#TH)/B 


D=2*CATNCQ1/SGR¢-Q1*Q1+1)>) 


POR SR—P7ie 70 ) STEP 07190 
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T5=TANCK/2) 
M=T3#T4/TS 
Mi=M*MO 

LPRINT " "3K 3M3M1 
NEXT K 

LPRINT " " 

aeeiNn " " 

aeeiNt * * 

LPRINT " " 


END 
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